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Abstract 

We present new near-infrared (IR) spectra (0.80-1.35/xm) of the pre-Main Se- 
quence source PV Cep taken during a monitoring program of eruptive variables 
we are conducting since some years. Simultaneous photometric and spectroscopic 
observations are systematically carried out during outburst and quiescence periods. 
By correlating extinction- free parameters, such as HI recombination lines and un- 
derlying continuum, it is possible to infer on the mechanism(s) responsible for their 
origin. Accretion and mass loss processes have a dominant role in determining the 
PV Cep irregular variability of both continuum and line emission. The potentialities 
of the observational modality are also discussed. 
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1 Introduction 



Young stellar objects (YSO's) present a certain degree of irregular variability 
which is typically attributed to intermittent accretion (or mass loss) phenom- 
ena. According to a widely accepted picture (Hartmann & Kenyon, 1985; Shu 
et al. 1994), accretion matter viscously migrates through the circumstellar disk 
toward its inner edge from where this matter is intermittently channeled along 
the magnetic field lines to the central star. The fall onto the stellar surface 
produces a shock that cools by emitting a hot continuum; moreover, as a con- 
sequence of the accretion event, strong winds (in some cases also collimated 
jets) emerge from the rotating star/disk system. Accretion-related variability 
may be due also to dust evaporation and condensation, a phenomenon po- 
tentially relevant in V1647 Ori (Reipurth & Aspin 2004, Acosta-Pulido et al. 
2007, Aspin 2011); V2492 Cyg (Hillenbrand et al. 2012); and PV Cep itself 
(Kun et al. 2011). Irregular variability is alternatively attributed to a variable 
extinction caused by dust clumps that occasionally intersect our line of sight 
(Skrutskie et al. 1996). Both mechanisms are plausible because all YSO's are 
expected to accrete, although with decreasing rates while the evolution goes 
on, and they are preferentially located in the densest parts of molecular clouds 
where gas and/or dust clumps are usual ingredients. 

The classical observational method for studying the modalities and the causes 
of the irregular variability consists of monitoring the continuum emission in 
different spectral bands to obtain the light-curves and the colors of the fluc- 
tuations. A further (and even better) technique is to exploit simultaneously 
a photometric and spectroscopic monitoring in the same spectral band: this 
represents a powerful tool of investigation since provides extinction-free pa- 
rameters that allow us to understand (i) the mechanism(s) responsible for 
the continuum and line emission; (ii) whether or not these mechanisms are 
(inter-)related; and (Hi) in affirmative case, to infer on the possible common 
mechanism. 

Remarkably, since few years we are conducting a photometric (Lorenzetti et 
al. 2006, 2007) and spectroscopic (Lorenzetti et al. 2009) monitoring in the 
near-IR (0.8-2.5 /xm, ZJHK bands) of a sample of low-mass pre-main sequence 
sources named EXor's (Herbig 1989; Reipurth & Aspin 2010). These are tar- 
gets very suitable for the aim described above, since they present large (about 
3-4 mag) outbursts, characterized by a short (weeks, months) rising and a 
longer (months, years) declining, superposed to long (years) quiescence peri- 
ods. In particular, the source PV Cep was extensively monitored (Lorenzetti et 
al. 2011; see also for a detailed description of this pre-Main Sequence object); 
its variability, firstly interpreted as due to extinction, has been later reconsid- 
ered as due to accretion and our aim here is to demonstrate (in the specific 
case of PV Cep) the effectiveness of the proposed approach. The novelty of 



2 



the present short note relies on both presenting new photometric and spec- 
troscopic data and, most importantly, considering new diagnostic parameters. 
Obviously, the approach adopted here for a single object (PV Cep) and for 
a relatively short monitoring period (about 3 years) can be ported in a more 
general context, once provided a meaningful monitoring of other sources. 

This paper is organized in the following way: Sect. 2 shortly refers to the used 
instrumentation; the relevant parameters derived from previous and new ob- 
servations are summarized in Sect. 3; our considerations are presented in Sect. 4, 
while final remarks are given in Sect. 5. 



2 Near-IR photometry and spectroscopy 

The following analysis is based on near-IR photometric and spectroscopic 
data. These latter were obtained at the 1.1m AZT-24 telescope located at 
Campo Imperatore (L'Aquila - Italy) equipped with the imager/spectrometer 
SWIRCAM (D'Alessio et al. 2000), which is based on a 256x256 HgCdTe 
PICNIC array. Photometry is performed with broad band filters J (1.25 /mi), 
H (1.65 /mi), and K (2.20 /mi). The total field of view is 4.4x4.4 arcmin 2 , 
which corresponds to a plate scale of 1.04 arcsec/pixel. All the observations 
were obtained by dithering the telescope around the pointed position. The 
raw imaging data were reduced by using standard procedures for bad pixel 
removal, flat fielding, and sky subtraction. 

Low resolution (1Z ~ 250) spectroscopy is obtained by means of two IR grisms 
Gbiue and G re d covering the ZJ (0.83 - 1.34 /im) and HK (1.44 - 2.35 /im) bands, 
respectively, in two subsequent exposures. The long slit is not orientable in 
position angle, and it samples a pre-defined portion of the focal plane, 2x260 
arcsec 2 in the east- west direction. 

Long-slit spectroscopy was carried out in the standard ABBA' mode with a 
total integration time ranging between 800 and 1200 sec. The spectral images 
were flat-fielded, sky-subtracted, and corrected for the optical distortion in 
both the spatial and spectral directions. Telluric features were removed by 
dividing the extracted spectra by that of a normalized telluric standard star, 
once corrected for its intrinsic spectral features. Wavelength calibration was 
derived from the OH lines present in the raw spectral images, while flux cali- 
bration was obtained from our photometric data carried out on the same night 
of the spectroscopy. 

During the monitoring period presented here, the near-IR spectra of PV Cep 
were taken in several occasions marked in Figure [T] with red vertical dotted 
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lines. Since this object is a target included in a photometric and spectroscopic 
monitoring program we are performing since few years, some of these spectra 
have been already discussed elsewhere (Lorenzetti et al. 2009). The oserva- 
tional novelty is represented by the 3 new spectra (and relative photometry) 
obtained in April 2010 (MJD 55310), July 2010 (MJD 55397), and August 
2010 (MJD 55416) (see Figure [T]). They will be commented, along with previ- 
ous data, in the following Sect. 4. 



3 Results 

The ZJ portion of the newly obtained (during 2010) near-IR spectra of PV 
Cep is plotted in Figure [2] along with our previous results. New data substan- 
tially confirm the existence of some continuum variability. HI recombination 
lines of the Paschen series are present: among them the Pa/3 appears as the 
most prominent in all spectra and on that line we will focus our attention. In 
Table [TJ the Pa/3 observed flux, its equivalent width (EW) and the value of 
the underling J band continuum (in mag) are given for each date. Given the 
spectral closeness between the J-band effective wavelength (1.25 /xm) and the 
Pa/3 transition (1.28 /an), the EW of this latter can be considered unaffected 
by the extinction. As such, EW is a suitable parameter to give indication on 
whether or not, and to what extent, the mechanisms responsible for the con- 
tinuum and line emission are related (see next Sect.). The same is obviously 
true also for other pairs of HI recombination lines/underlying continuum: in- 
deed we did the same analysis also for Br7/K band combination obtaining the 
same results not reported here to avoid any unuseful redundance. 



4 Analysis and discussion 

In the upper panel of Figure [3] an evident correlation (regression coefficient = 
0.93) is recognizable between the Pa/3 line flux and the continuum in the J 
band. Such a correlation suggests, as expected, that emission lines are gener- 
ated by some mechanism strongly related to the central source. According to 
the lower panel of the same figure an ant i- correlation seems to exist between 
EW and J continuum, in the sense that decreasing EW values correspond to 
an increasing continuum brightness. Two data points, those corresponding to 
J about 12 mag and marked with open symbols, seem not to follow the anti- 
correlation trend, but this apparent contraddiction will be discussed later. 
Neglecting these two points, a regression coefficient of 0.91 is obtained. Now 
we remark the anti-correlation indicates that, although line emission and J 
continuum are related, the latter presents an increase larger (i.e. faster) than 
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Table 1 

Spectroscopic and photometric parameters of PV Cep. 



Date 


MJD" 


Pa/3 Flux 


EW 


J band 






(10~ 13 ergs _1 cm~ 2 ) 


(A) 


(mag) 


May 13, 2007 


54241 


3.8 ± 0.2 


-48 ± 2.4 


11.52 6 


Oct 15, 2007 


54389 


3.7 ± 0.2 


-40 ± 2.0 


11.53 


Oct 27, 2007 


54401 


3.5 ± 0.2 


-28 ± 1.7 


11.14 


Apr 01, 2008 


54558 


5.3 ± 0.2 


-16 ± 0.6 


10.09 


Apr 12, 2008 


54569 


6.1 ± 0.2 


-20 ± 0.6 


10.14 


Jun 10, 2008 


54627 


1.4 ± 0.2 


-22 ± 3.0 


11.96 


Jun 18, 2008 


54635 


1.2 ± 0.2 


-19 ± 3.0 


12.34 


Apr 24, 2010 


55310 


5.5 ± 0.2 


-29 ± 1.3 


10.63 


Jul 20, 2010 


55397 


5.1 ± 0.2 


-25 ± 1.6 


10.58 


Aug 08, 2010 


55416 


5.9 ± 0.2 


-31 ± 1.8 


10.63 



a MJD = Modified Julian Day 

b Errors on J magnitudes are typically 0.02-0.03 mag 

the former. Indeed, the J-band flux is composed of the direct radiation from the 
hot spot and an indirect brightening due to the radiation from the innermost 
regions of the circumstellar disk. Since PV Cep is embedded in an extremely 
massive disk (Hamidouche 2010), seen at large inclination, this effect has to 
be accounted for and likely represents the additional contribution that makes 
the continuum variation larger than the line emission. Similar anticorrelation 
results are given by Cohen et al. (1981), Magakian & Movsessian (2001), and 
Acosta-Pulido et al. (2007). 

This occurrence tends to rule out the variable extinction as cause of the ob- 
served variability that affects to the same extent both lines and continuum, 
since in that constant value of the EW should be expected for any con- 

tinuum fluctuation. Even more a selective obscuration can be also ruled out. 
In fact, it should be located at the dust condensation zone, very close to the 
star and, as such, it should obscurate the stellar photosphere more than the 
accretion columns and the wind regions. On the contrary, accretion (or mass 
loss) processes are compatible with the observed behaviour, since they respond 
with different laws (Muzerolle et al. 1998), and maybe with different lags, to 
the continuum variations. These latter, once ruled out any extinction varia- 
tion are unanimously attributed to hot spots on the stellar surface due to the 
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shocks that originate at the base to the columns of the infalling matter. Being 
the line fluxes well correlated with the increasing continuum, they will depend 
on accretion processes: such a dependence may be direct if lines originate in 
the accretion columns, or indirect if lines originate in a stellar wind whose 
mass loss rate (M win d) is related to the mass accretion rate (M acc ) by invoking 
the rough proportion M win d/M acc < 0.1 (Shu et al. 2000; Konigl & Pudritz 
2000). Therefore we can assume that the total HI recombination line flux is the 
sum of different contributions: from accretion, from wind and, possibly, from 
other components (e.g. an HII circumstellar region). During strong accretion 
events, that increase the source continuum emission, only one contribution to 
the line emission dominates (correlation depicted in Figure [3]), but when the 
impulsive accretion stops, namely when the continuum drops below a given 
threshold (in the present case J mag > 12), line emission strength comes 
back to its original (lower) value, practically independent on the continuum 
value, and the EW value suddenly changes (see open symbols in Figure [3J. 
In principle, alternative interpretations could exist: the proposed one is well 
compatible with the observational data and reconciles the two apparently un- 
consistent data points. However, to firmly confirm the presented scenario we 
still needs additional data resulting from the photometric and spectroscopic 
monitoring of relatively strong (more then 1-2 mag) variations of PV Cep or 
other targets. Indeed, one aspect that would deserve a deeper investigation 
is the following: by looking at Figures [4] and [TJ we can realize that the EW 
value reached immediately after the accretion event (that indicated by open 
symbols) is very similar to that corresponing to observations taken 100 days 
before, at an outburst peak. According to the proposed scenario, any EW 
value is associable to a continuum value lesser than a threshold, hence the 
observed EW resemblance is remarkable, but, plausibly, fortuitous. 



5 Final remarks 



The case presented here for a single young source (PV Cep) aims at fos- 
tering systematic photometric and spectroscopic monitoring intended as an 
observational method of investigation on how the mechanisms responsible for 
continuum and line emission are (inter-)related. Unfortunately, the proposed 
technique is not widely exploited since a frequent (weekly) monitoring (on 
years time-scale) is not easily doable both on large telescopes affected by a 
heavy oversubscription and on smaller telescopes often not equipped with the 
needed instrumentation. It would be very fruitful if small telescopes (l-2m 
class), but properly equipped, were dedicated to similar projects even for a 
small fraction (< 10%) of their observing time. 
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Fig. 1. PV Cep near-IR light curves vs. MJD (Modified Julian Date). The errors of 
data points are comparable to or lesser than 0.03 mag. Vertical red lines correspond 
to those dates when, beside JHK photometry, also a near-IR spectrum was taken. 
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Fig. 2. The newly obtained near-IR spectra of PV Cep are depicted in colours (red, 
green, magenta). For comparison purposes, the already presented spectra are also 
plotted (in black). 
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Fig. 3. Upper panel. Distribution of the Pa/3 fluxes extracted from the different 
spectra of PV Cep as a function of its J magnitude. Lower panel. Same as before, 
but in terms of Equivalent Widths; as usual, EW have negative values because line 
are in emission. Open symbols refer to the spectra taken just after a rapid fading 
of the source (see text and Figure H| . 
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Fig. 4. EW derived from Pa/3 spectroscopy as a function of time. As in Figure [3] 
the same couple of data points is represented with open symbols. 
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